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How do we uniquely identify a quantum phase, given its ground state wave-function? This is a 
key question for many body theory especially when we consider phases like topological insulators, 
that share the same symmetry but differ at the level of topology. The entanglement spectrum has 
been proposed as a ground state property that captures characteristic edge excitations. Here we 
study the entanglement spectrum for topological band insulators. We first show that insulators with 
topological surface states will necessarily also have protected modes in the entanglement spectrum. 
Surprisingly, however, the converse is not true. Protected entanglement modes can also appear 
for insulators without physical surface states, in which case they capture a more elusive property. 
This is illustrated by considering insulators with only inversion symmetry. Inversion is shown 
to act in an unusual way, as an antiunitary operator, on the entanglement spectrum, leading to 
this protection. The entanglement degeneracies indicate a variety of different phases in inversion 
symmetric insulators, and these phases are argued to be robust to the introduction of interactions. 



I. INTRODUCTION 

According to the Landau paradigm, phases of mat- 
ter are classified by their pattern of symmetry break- 
ing. While this accounts well for a number of experimen- 
tally observed phenomena, several exceptions have also 
cropped up. These include the fractional quantum Hall 
states and models of gapped spin liquids, which are char- 
acterized by a topological order: for these phases, the 
ground state degeneracy depends on the spatial topology 
of the sample. Perhaps the simplest phases for which 
a topological distinction is present are topological insu- 
lators and superconductors, composed of noninteracting 
particles. Experimental manifestations include the inte- 
ger quantum Hall effect, spin-orbit induced Z2 topologi- 
cal insulators in twcP^ and three dimensions,^ ^ as well 
as superfluid He^ — B. Here, although there is a unique 
ground state and no symmetry difference from the triv- 
ial state, the topological distinction can manifest itself 
in different ways. These include a quantized response 
functiorP (such as the Hall effect), or protected surface 
states. Given the subtleties associated with identifying 
topological phases, a new tool in this respect would be 
welcome. 

Recently, it has been shown that studying the entan- 
glement spectrum is a promising direction to identifying 
topological phases. Given the ground state wavefunction, 
and a partition of the system into a left and right half, 
one can perform a Schmidt decomposition: 

\G) = Y.^e-^\aL)\aR) (1) 

Measuring the right half of the system shows it to be 
in state \aR) (called a Schmidt state) with probability 
e~^«/Z. The quantities comprise the entanglement 
spectrum, and are somewhat like "energies," character- 
izing how unlikely a given fluctuation is to occur. (Note 
that we use the superscript "e" when describing the en- 
tanglement.) The entropy associated with this probabil- 



ity distribution is the entanglement entropy, which cap- 
tures the same features of a phase as ground state degen- 
eracy (or topological order).^ However, the entire spec- 
trum can be expected to contain more information, al- 
lowing one to capture a wider class of distinctions. Even 
though entanglement related quantities are hard to mea- 
sure directly, they can be conceptually very useful in 
identifying phases. 

It has been shown for fractional quantum Hall states 
that although their bulk is gapped, the entanglement 
spectrum contains information about the edge modes. 
Numerical studieJ^^ reveal that the largest entanglement 
eigenvalues (smallest E^) mimic the low-energy spectrum 
of surface modes. A relationship of this type is very in- 
teresting because the entanglement spectrum can be cal- 
culated entirely from the ground state wave function. In 
contrast, the edge modes are excited states of the Hamil- 
tonian in a sample with boundaries. In a classical sys- 
tem, one would not be able to look at the ground state 
and determine its dynamics. But a quantum mechani- 
cal system, even in its ground state, has some zero-point 
motion that can give one a sense of its excited state prop- 
erties. A direct classification of phases from their ground 
state wavefunctions is highly desirable and could be ap- 
plied to cases where a candidate ground state (such as a 
Gutzwiller projected state) is available. 

Here we study the entanglement spectrum in the con- 
text of the simplest class of topological phases - noninter- 
acting topological band insulators. We show that when- 
ever topological surface states are present, the entangle- 
ment spectrum, as a function of the momentum parallel 
to the cut, also has protected entanglement modes. Es- 
sentially, the entanglement spectrum is known to result 
from diagonalizing an operator, and this operator can be 
viewed as a band insulator Hamiltonian which retains the 
topology of the physical insulator. Explicit examples of 
the entanglement spectrum are worked out, including an 
analytical calculation for Landau levels (similar to the 
work of Ref.TT") and a numerical study of three dimen- 
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sional topological insulators. Similar reasoning can be 
applied to topological superconductors. Such a relation 
has also been noticed in isolated examples, like the quan- 
tum Half^ and the p^^+zp^y superconductoi^^ edge. Inde- 
pendent studies^^ have appeared, pointing to a similar 
connection. 

Next, we consider the converse - i.e. do protected 
modes in the entanglement spectrum necessarily imply 
protected surface excitations? We show explicitly that 
this is not true, by studying an example of a insulator 
with only inversion symmetry. Here no protected sur- 
face excitations exist, but the entanglement spectrum 
features protected states, which points to a more sub- 
tle distinction between phases. The entanglement spec- 
trum remains gapless because inversion symmetry is re- 
tained when dividing the system for the purposes of cal- 
culating entanglement entropy, but is implemented in a 
strange way, as a particle-hole symmetry of the single 
particle entanglement states. This leads to protected 
modes. Thus studying systems with inversion symmetry 
(or another symmetry that maps the left and right halves 
to one another) is a good way to break the close corre- 
spondence between edge states and entanglement, a point 
noticed earlier for one-dimensional interacting states in 
Ref. I15|16[ When this correspondence breaks down, the 
protected entanglement degeneracies still indicate a dis- 
tinct phase, although there are no physical surface states. 
Understanding the consequences of this hidden degen- 
eracy should be interesting. Although disorder breaks 
inversion symmetry, and is hence not normally studied 
while classifying topological insulators j -*^^ * "*^^ * clean phys- 
ical systems with inversion symmetric bulk states can 
certainly be realized. 

Finally, we recast the theory given here in a form that 
is suitable for studying interacting systems: the action of 
inversion (/) is defined on the many-body Schmidt states. 
It is shown to act as an anti-unitary operator /M, which 
must satisfy (/M)^ = +1 or (/M)^ = -1. The latter cor- 
responds to the topologically nontrivial case, and leads 
to a two fold degeneracy of all states in the entanglement 
spectrum. This degeneracy remains on introducing weak 
interactions. Part of the topological distinction implied 
by this observation is contained in the quantized mag- 
netoelectric polarizability of insulators, which remains 
quantized in inversion-symmetric systems. We also dis- 
cuss finer distinctions between inversion-symmetric insu- 
lators, but their physical implications remain to be iden- 
tified in future research. 



II. TOPOLOGICAL BAND INSULATORS 

A band insulator is described by the single particle 
Hamiltonian: 

if = ^e^(k)dy^k (2) 

7,k 



where k is crystal momentum, and 7 = 1,...,A/' is a 
band index, of which n bands are filled i.e. e^(k) < 
if 1 < 7 < n and e^(k) > if n < 7 < TV. The band 
wave-functions are: 

4k = E'^Tk(r)V'^(r), (3) 

r 

where we have suppressed spin and orbital indices (These 
variables can be included in r with the spatial coordi- 
nates). 

Topological phases of band insulators are defined by 
their band topology. Two insulators differ topologically 
if they cannot be connected by smooth changes of the 
Hamiltonian while the band gap remains finite. It is 
possible to define a topological invariant to distinguish 
between different phases that depends only on the wave- 
functions of the filled bands. ^ 

Given a general band insulator, ([2|, a topologically 
equivalent insulator can be constructed by setting the 
energy of all occupied bands to be equal and negative, 
and all unoccupied bands to be equal and positive: e.g. 

CF^ik) = -iif7<n (4) 

cfjO^) = +^if7>n. (5) 

Let us call the corresponding operator Q. Note, this is 
related to the projection operator P on the filled bands 
via Q = ^1 — P. This "fiat band" limit has been found to 
be useful in the classification of topological insulators 
and will be used here as well. An important connec- 
tion is with the correlation function C(r,r') = (V^^'^r), 
evaluated in the ground state. Now P = (7, i.e. the 
correlation function of the band insulator, viewed as a 
matrix, is simply the projection operator onto the filled 
bands. Hence 

Q = E I'^rA - E ^(r- ^r' . (6) 

r 

describes the fiat-band insulator. Since correlations in a 
gapped state like a band insulator fall off exponentially 
rapidly with separation, the effective "flat band" Hamil- 
tonian (|6| has essentially short ranged matrix elements, 
as for a physical operator. Hence, it can be viewed also 
as a bona fide Hamiltonian, with the same band topology 
as the starting Hamiltonian in Eqn. [2] 

Most band insulators with nontrivial band topology 
are characterized by edge states. Consider representing 
the Hamiltonian of a band insulator in real space: 

H = ^^(r,r')V^(r)V(r') (7) 
n{r,r') = J2e^(k)<j>;^ir)Mr') (8) 

7k 

A boundary along the plane x = 0, with the physical 
system to the right {x > 0), is obtained by truncating 
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FIG. 1: A cut that divides the physical system into two halves. 



If both the system and the cut preserve translation 
symmetry, one can simultaneously diagonalize the den- 
sity matrix and the translation operators parallel to the 
cut Ty. The entanglement spectrum is then obtained 
as a function of momentum along the cut: £^^(Kx). 
This object can capture fairly detailed properties of the 
ground state wavefunction that are discussed in what fol- 
lows. 

For the case of single particle Hamiltonians, such as 
Eqn. [2j Ingo Peschel has shown how to determine the en- 
tanglement Hamiltonian and its eigenvalues.^^ The first 
step is to realize that the entanglement Hamiltonian must 
be quadratic in the Fermionic operators: 



the Hamiltonian: 

H^(r,rO = Oifx<0orx'<0 (9) 
H^(r,r') = H(r,r') otherwise (10) 

The spectrum of states localized deep in the bulk is 
unaffected by the cut, due to the presence of a finite 
gap. However, as one approaches the edge, it is possible 
that states appear within the bulk gap. Most topolog- 
ically nontrivial insulators have "protected" edge states 
throughout the bulk gap; they cannot be removed by 
smooth changes of the Hamiltonian. Thus, if the Hamil- 
tonian ([2| has protected edge modes, so must its "flat 
band" version (|6|), when restricted to a half space. 



III. THE ENTANGLEMENT SPECTRUM 

Now we will review what the entanglement spectrum 
is, and show that the edge modes of the flat band Hamil- 
tonian are in fact edge modes of the entanglement spec- 
trum. 

Consider the ground state wavefunction of a general 
Hamiltonian H, "^ciia^}, {o^f})^ where {af } ({af }) 
specifies the configuration of the system to the right (left) 
of an imaginary cut. Tracing out the degrees of freedom 
on the left, one obtains the density matrix: 

PR{{af}, {«f }) = E *G({af }, {af })*2;({af }, {af }). 

(11) 

For a normalized wavefunction, Tt[pr] = 1, and the 
eigenvalues of pR lie in the range [0, 1]. Thus we can 
write: 



m 



(13) 



PR = Z-'^e-""^; Z = Tre-""^. 



(12) 



We will call Hf^ the entanglement "Hamiltonian". The 
eigenvalues of the density matrix are denoted by pa = 
Z~^e~^«, where the set is termed the "entanglement 
spectrum". Note, these "energies" have nothing to do 
with the physical energies of the system - they refer en- 
tirely to properties of the ground state. (Note that Z is 
introduced in order to allow us to set the zero of "energy" 
conveniently.) 



This can be seen, and the single particle "entanglement 
Hamiltonian" determined, by noting that any corre- 
lation function of operators that only involve degrees of 
freedom on the right is unaffected by tracing out the left 
half of the system. For example, consider the correlation 
function C(r,r') = (?/^J/?/^r), in the ground state, before 
the system is cut. Now, if both coordinates r, belong 
to the right half of the system, then they are unaffected 
on tracing over the left half. Therefore, the same result 
should be obtained by using the density matrix for the 
right half. 



(14) 



where Cr is the correlation function restricted to this 
half, namely Cr{y^y') = C(r, r') if both coordinates are 
on the right, and zero otherwise. A similar argument 
applies to multiparticle correlations. Thus, since Wick's 
theorem is obeyed by correlation functions in the orig- 
inal uncut system, it will continue to hold for just the 
right half of the system in the mixed state p^. Conse- 
quently, p^ must be Gaussian, i.e. the exponential of a 
quadratic form of Fermion operators. Requiring also con- 
servation of particle number leads to the general form in 
Eqn. 13 An expression for may be obtained by re- 



turning to two-point correlations, and requiring that the 
"Boltzmann" distribution, Eq. (12), gives the expres- 
sion Eq. (14). One expands the Fermion operators in 
a basis which diagonalizes H^; in this basis, C becomes 
a diagonal matrix recording the mean occupation num- 
bers of the states, given by the Fermi distribution. Thus, 
CR = [l^e^'']-\oi 



log 



1-Cr 
Cr 



(15) 



Equivalently, the operator Qr = ^1 — Cr^ satisfies the 
relation: Qr = ^tanh[^H^]. Thus if Pi are the eigenval- 
ues of Cr^ which are readily seen to satisfy < < 1, 
then the eigenvalues of the 'single particle' entangle- 
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ment Hamiltonian are obtained from Eqn 15 



el = loe 



Pi 



(16) 



The many body entanglement spectrum is obtained 
by deciding the occupancy of the single Fermion modes, 
so that mode i has nf =0, 1. Then = e^. The 

density matrix eigenstate with the largest eigenvalue is 
like the Fermi sea of a physical system; it corresponds to 
the lowest "energy" E^ obtained by filling all ef < 0. 

For a translationally symmetric cut the entanglement 
energies can be resolved as a function of the total trans- 
verse momentum (following Ref. 10). These energies 
can again be built from the single particle "energies" of 
the occupied states, e^(k^), which are functions of the 
momentum parallel to the cut. The total momentum 
Kx is also just the sum of the momenta of the filled sin- 
gle particle states. Thus, the fundamental object we will 
focus on calculating is the single particle "spectrum" as 
a function of transverse momentum e^(k±). 



IV. CONNECTING THE ENTANGLEMENT 
SPECTRUM TO EDGE MODES 

The previous discussion clarifies the precise connec- 
tion between the entanglement spectrum and the edge 
spectrum of an insulator. The single particle entangle- 
ment spectrum ef(kx) is related to the spectrum ob- 
tained when a flat band version of a physical Hamiltonian 
is studied in the presence of an edge. The eigenvalues of 
the flat band Hamiltonian Q, which are ^ — Pi(k±)^ are 
related to the entanglement spectrum ef(kx) via Eqn. 
(IT6|), or 



^tanh[ie,^(k^) 



(17) 



Most states have eigenvalues ±oo, because of the 

nonlinear relationship; the eigenvalues of C in the bulk 
are 1 and 0. Very negative eigenvalues represent bulk 
states that are occupied, and large positive eigenvalues 
represent unoccupied stated. 

In cases where topologically protected surface states 
of the physical Hamiltonian are expected, the fiat band 
deformation which is topologically equivalent is also ex- 
pected to have surface states^ filling the whole gap be- 
tween the bulk states at +oo and — oo. 

Note however, the converse is not necessarily true. For 
inversion symmetry (see below), the edge of the flat band 
Hamiltonian has structure which signifies protected fea- 
tures of a phase, although this structure is not present 
in the physical surface states of a generic Hamiltonian in 
the same topological class. 

We will now illustrate our reasoning with various ex- 
amples. 




FIG. 2: Entanglement edge states for the v — 2 integer quan- 
tum Hall state. 



V. EXAMPLES 

Let us summarize the procedure for obtaining the en- 
tanglement spectrum. The following prescription takes 
advantage of k^ conservation to reduce the problem to a 
one-dimensional problem: (i) obtain the correlation func- 
tion restricted to the right half, by summing over occu- 
pied bands. If the wavefunction of the filled orbit als at 
momentum k = {kx^ k^) is (/)ike*^ '^, where r is the unit 
cell position, and a refers to other indices such as sublat- 
tice and spin, then: 

C^^{x,a;x',a') = ^ e''===(---'),^*k(«')'/'ik(a). 



(ii) Find Cs eigenvalues ^^(k^), by solving the eigenvalue 
equation 

Y^C^S^, a-x\ a')fg^{x\ a')=p^{k^)fgjx, a) 



(iii) Then the single particle entanglement spectrum 
ef (k^) can be read off from equation (17) and the eigen- 
functions are fi\^j_{x)e^^-^'^-^ . We will graph our results 
for ^ —p^(k±) instead of e^(k^), since they are monoton- 
ically related. 



A. Integer Quantum Hall: 

The first example is based on the work of Ref. TP on 
the integer quantum Hall effect. This is an especially 
clear illustration of the general relationship between the 
correlation function and the entanglement Hamiltonian. 
Consider electrons confined to a plane and in a strong 
magnetic field; assume there is no crystal potential and 
suppose that exactly u Landau levels are fully occupied. 
We show that the entanglement spectrum contains the 
expected u chiral edge modes. Using the Landau gauge, 
Ay = Bx, and labeling the states by the momentum 
along the cut k± = ky^ the restricted correlation func- 
tion is given by: C^^{x, x') = ELi C/c, (^O^n/c, (2:) 
where the normalized eigenstates are (j)nky{x) = Hn(X — 
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Ky)e Jiere X = x/l; Ky = kyl are dimension- 

less variables scaled by the magnetic length / = ^/h/eB^ 
and Hn are appropriately normalized Hermite polynomi- 
als. 

Now, an eigenfunction of must be constructed 
from linear combinations of the Landau level states = 
Ylm^i^m^m (whcrc ky has been suppressed). The coef- 
ficients are easily seen to be eigenvectors of the u x u 
matrix Fnm{ky) = dx (j)lf^^{x)(j)mky{x). For = I, 
this is a number, which is the eigenvalue itself. Thus 

l/2-po{ky) = -^evi{Ky), 

which interpolates between —1/2 (when ky = — oo) and 
-hl/2 (when ky = +00). This is the single chiral mode, 
as expected. Similarly, the u = 2 case can be solved 
analytically, the result for the two eigenvalues is now 

leading to two modes as shown in the figure. These two 
modes are like the dispersion of two chiral particles in 
a one-dimensional system. In general, the entanglement 
Hamiltonian of a gapped system behaves like the Hamil- 
tonian of a system with one dimension fewer. There may 
be infinitely many bands besides those near = 0, but 
their energies rapidly approach ±00. 




FIG. 3: Brillouin zone of the (1,1,1) surface of the diamond 
lattice. The labelled points are time reversal invariant mo- 
menta. 

r — To ^ To — r, about a bond center tq). Without the 
second-neighbor hopping, the model's bulk and surface 
spectra display the non-generic feature that all energy 
levels come in pairs at each momentum. 
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B. Three Dimensional Topological Insulator with 
Time Reversal Symmetry: 

We now calculate the entanglement spectrum for a 
three dimensional topological band insulator, an insu- 
lator with surface modes protected by time reversal sym- 
metry. By explicit calculation we confirm the expectation 
that these modes are also captured by the entanglement 
spectrum, both in the case of strong and weak topological 
insulators. 

We consider the model of a topological insulator on 
the diamond lattice introduced in Ref.4, given by the 
Hamiltonian: 



Sit SO 



'CkJ^lS) 



{{ik)) 



where the first term contains the four nearest neighbor 
hopping elements which are taken as t for three of the 
bonds, and t-\-5t for the fourth bond oriented along the (- 
1,1,-1) direction. To access a generic Hamiltonian we also 
include regular second-neighbor hopping ^2 (not shown in 
Eq. 18) in all directions. The spin orbit interaction tso 
appears in the second term, inducing hopping between 
second neighbor sites. In this term, d-^^, d-^^ are the two 
nearest neighbor bond vectors leading from site i to /c, 
and a are the spin Pauli matrices. This model respects 
time reversal symmetry, and also inversion symmetry (i.e. 



M1 



M3 



(b) 



FIG. 4: Strong topological insulator, a) The dispersion of 
the physical surface states. Note that the node is not at zero 
energy, since particle-hole symmetry is broken, b) The entan- 
glement spectrum, | — p, which has the node at zero e^, in 
spite of particle-hole symmetry being broken. 



We now compare the physical surface states that ap- 
pear in this model and the entanglement spectrum of 
a cut with the same orientation. We choose a surface 
normal to the (1,1,1) direction so as to cut bonds with 
strength t. One obtains a strong topological insulator, 
with an odd number of surface Dirac nodes if 6t > 
(assuming t > 0). The surface states computed for a sin- 
gle surface are shown in Figure Qa), and a single Dirac 
node characteristic of the strong topological insulator is 
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obtained. The node is centered at the time reversal in- 
variant M2 point of the surface Brihouin Zone (BZ) (see 
Figure [3|. This is selected by the direction of the strong 
bond t + St. Note that the surface spectrum is not sym- 
metric between positive and negative energies, as it would 
be if the model were particle-hole symmetric. 

Below, in Figure |4jb), we display the two surface 
eigenvalues of the single particle entanglement spectrum 

(k^), obtained from the ground state of the Hamilto- 
nian in Eqn. [l8]by dividing the system into two halves. 
The plane that divides the system is oriented in the same 
way as the physical boundary previously discussed. It is 
more convenient to display ^tanh(^^^^^) which is re- 
lated to the eigenvalues Pi(k±) of the correlation func- 



tion Cr via ^ tanh( 
ergy", e^^{k± 



2 J — 2 PiO^-^)' Note zero "en- 
0, also corresponds to the zero of ^ 



-Pi 



Clearly, we see that the entanglement spectrum also dis- 
plays the same characteristics as the protected surface 
modes. There is a single Dirac cone which is centered 
at the same point in the surface BZ as the surface state 
itself. The "energy" of this Dirac cone is curiously zero 
(understanding this is our next job). 

Similarly, if we choose < 0, we obtain a weak topo- 
logical insulator, whose physical surface states and entan- 
glement spectrum are compared in Figure |5) for the same 
surface as above. Again, the Dirac nodes of the physical 
surface states, and their location in the Brillouin zone, 
which are fixed by band topology, are captured by the 
entanglement spectrum. 

We now briefiy describe the computation that gives us 
the entanglement spectrum in the figures. Essentially, 
we follow the three step procedure outlined in Section 
\V\ where is a momentum on the surface BZ. Instead 
of dealing with a pair of semi-infinite systems L, it is 
more convenient to divide the system into three parts, a 
central piece, whose reduced density matrix is computed, 
and one part to the right and one to the left, which one 
traces over. This gives us two entanglement surfaces, 
but since they are well separated by a gapped region, 
the spectra are essentially independent, and we only plot 
those whose wavefunctions are concentrated near one of 
the cuts. For the numerics we have taken the separation 
to be 80 unit cells thick, and we have taken the following 
parameters in the Hamiltonian Eqn. (18): t = 1.0, tso = 
0.125, St = 1 (-0.9), t2 = 0.1 for the strong (weak) 
TI. The bulk energy bands as well as the entanglement 
spectrum eigenvalues near 1/2— p = ±1/2 are not shown 
for clarity. 

This has demonstrated that when a protected physical 
surface state is present, this is refiected in the entangle- 
ment spectrum. That the Dirac nodes in the entangle- 
ment spectrum are centered at zero, in contrast to the 
physical surface spectrum in the same model, is an im- 
portant observation. Moreover, the entire entanglement 
spectrum remains symmetric under — e^, k^ 

— kx- We will see that this is a consequence of inver- 
sion symmetry, and the Dirac node in the entanglement 
spectrum remains secure even when time reversal invari- 




(a) 




(b) 



FIG. 5: Weak topological insulator, a) The dispersion of the 
physical surface states, b) The entanglement spectrum. 



ance is broken. This is shown in Figure |6] The physical 
spectrum of the surface states Fig. [6^ is gapped because 
of a uniform Zeeman field Ht = h^^c\(7zCi applied to 
all the sites, with h = 0.9. However, the entanglement 
spectrum Dirac node remains intact, despite the loss of 
time reversal symmetry in Fig|6]3. When both time re- 
versal symmetry and inversion symmetry are destroyed, 
the latter via a staggered potential Hj = V ^-{—lyc^ Ci 
(with V = 0.1), then the entanglement spectrum is also 
gapped as in Figure [TJd. 



VI. TOPOLOGICAL INSULATORS WITH ONLY 
INVERSION SYMMETRY 

We now discuss the origin of Fig. [6j where a pertur- 
bation breaking time-reversal symmetry has been added. 
The entanglement spectrum still has a Dirac node, which 
is in fact at zero entanglement "energy", although the 
surface states are gapped. Therefore, while a system with 
topologically protected surface states also has topologi- 
cally protected entanglement states, the converse may 
not be true. The Dirac node is a feature in the entangle- 
ment spectrum of inversion symmetric insulators, which 
defines a distinct phase, although there are no topolog- 
ical surface states for the physical Hamiltonian generi- 
cally. The entanglement spectrum is an especially useful 
tool for identifying insulators of this type. 



r M1 M2 M3 r 

(a) 



r M1 M2 MS r 

(a) 




(b) 



FIG. 6: Inversion symmetric insulators with broken time re- 
versal symmetry. Note the physical surface states (a) are 
gapped, but the Dirac node remains in the entanglement spec- 
trum (b). 



The reason the entanglement spectrum works for iden- 
tifying phases with inversion symmetry, while physical 
surface states fail, is as follows. A physical boundary 
necessarily breaks inversion, since inversion maps the in- 
sulator to its exterior. But now consider a cut that passes 
through a center of inversion. The inversion X through 
this point relates the left and right sides of an entangle- 
ment cut. This inversion symmetry leads to a symme- 
try of the entanglement Hamiltonian, but with the twist 
that it is realized as a sort of particle-hole symmetry 
for the two dimensional entanglement spectrum (see also 
Refs.HSnnj). 

Let us suppose that the inversion symmetric cutting 
plane does not exactly coincide with any orbitals(so that 
it divides the states into separate parts)^. We call the 
inversion transformation through the center of inversion 
on the plane X. In this section, we will describe how 
this symmetry constrains the entanglement states, then 
use this to prove the stability of the entanglement Dirac 
node. Finally we discuss what physical difference could 
be embodied in this topologically distinct entanglement 
spectrum. The remaining sections will discuss interacting 
insulators. 

To see why inversion symmetry acts as a particle-hole 
symmetry we will give an additional interpretation of the 
entanglement eigenstates f^^. The entanglement spec- 
trum can be understood most intuitively by showing that 




(b) 



FIG. 7: Without inversion symmetry. Both the physical sur- 
face state (a) and the entanglement spectrum (b) are gapped. 



a set of wave functions for the full systenP^I^ 



= "v/^^/^ fx) 



v/i-p.(k^)/^k,(^) 



(19) 



can be constructed, satisfying two properties: first, the 
/^'s and /^'s form an orthonormal family of wave- 
functions in the two halves. Second, filling all the states 
^ik^ gives the band insulator; i.e. they can be thought 
of as being obtained from the filled band wavefunctions 
(Eq. [s]) by a unitary transformation. The /^'s and the 
Pi's turn out to be the same as before. (See the ap- 
pendix.) The eigenvalues ^^(kx) have a simple interpre- 
tation in light of Eq. ([l9|. Each state F^k^ is certainly 
occupied by an electron in the band insulator. The eigen- 
value Pi(k^) represents the probability that the electron 
in this state is found on the right half. For example, 
when ^ 1, corresponding to — oo, the electron 

in Fi is almost certainly on the right half. This fits with 
the previous definition of the /^'s because is far be- 
low the "Fermi energy," so the z*^ state is occupied in 
the "Fermi sea" of the entanglement Hamiltonian. Eqn. 
(19) leads to a very intuitive way of understanding the 



entanglement spectrum, as reviewed in appendix A. 

Because of the inversion symmetry, there is a map Xs 
on the single-particle states that takes a state to 
another state ffl_\^^- This state flips the sign of the en- 



8 



tanglement energy, e|(— k^) = —e^(k±). The mapping 
arises directly from inversion symmetry acting on the F's 
(see Fig. [S^i), which produces from any occupied state Fi 
with momentum a second state, with the opposite 
momentum: 



Expanding this equation gives 



(20) 



Since inversion maps the left-hand side to the right-hand 
side, these equations imply that k^) = 1 — Pi(k±) 
and that f-^ is obtained from a state on the other side of 
the partition, namely f^_^^{x) = f^^{—x). Using the 
relation between p and e, it follows that a mode with 
"energy" and momentum k^, is mapped by inversion 
to one with — and — k^. A more algebraic proof of this 
result is in Appendix [A) 

Having established the action of inversion, we now turn 
to the stability of Dirac modes in the entanglement spec- 
trum, in the presence of inversion. Suppose the insula- 
tor is obtained from a time-reversal symmetric topolog- 
ical insulator by applying a time reversal breaking per- 
turbation. (We plan to derive the basic prop erties of 
general inversion symmetric insulators later A Dirac 
mode in the entanglement spectrum (or physical surface 
spectrum) occurs at a time reversal invariant momentum 
(or TRIM) n^. These momenta are half of a recipro- 
cal lattice vector, i.e., k±_ G {F, Ml, M2, M3}, so that 
= —K± (modulo the reciprocal lattice). 

First consider the modes in the unperturbed crystal, 
that is both time-reversal and inversion symmetric. At 
the tip of the cone of a Dirac mode, there are two de- 
generate states /^^, /^^. These must have = 0, 
because otherwise applying Is would produce a second 
Dirac point at — e^. 

Now the symmetry Xs ensures that these Dirac modes 
have to remain at zero "energy" when time-reversal sym- 
metry is broken by a small perturbation. First, the 
two states are transformed into themselves under X^, 
so appropriate linear combinations of them are inver- 
sion eigenstates. In fact, their inversion parities must be 
the same. (The states formed a Kramers doublet under 
time-reversal symmetry before it was broken.) There- 
fore the states f^^_^ and cannot evolve into a pair 
f^,K±^ f-,K± with energies ±e^, without a discontinu- 
ous jump in the inversion parities. (Inversion symme- 
try interchanges these new states, so their combinations 
71 (/+,/^^ ^f-,K,J have opposite parities.) 

Note, in general, the spectrum consists of equal num- 
bers of states that are even under inversion and odd under 
inversion. All the pairs of states at nonzero entanglement 
energy can be combined into pairs of states of opposite 
parity, since /e,k^ ^ /-e,-k^7 and one can take sym- 
metric and antisymmetric combinations of these. The 
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FIG. 8: Inversion symmetry of entanglement, (a) The trans- 
formation Ts of single-body modes. The image of a state 
under Xs is not defined on the merits of itself (like ordinary 
reflection), but instead depends on the state of the system. 



The extended wave function (Fi — y/piff^ + \/l — p%ft } is m- 
verted, and Xsf^ is (up to normalization) the right-hand part 
of this, which can look completely different from the original 
(b) The transformation MX of many-body states in a 
non-interacting system. Each term in the Schmidt decompo- 
sition is obtained by placing electrons on either the right or 
left of the Fi^s indicated by the horizontal lines. Inversion is 
applied to the system as a whole. Focusing just on the half 
right-hand side of the system, one finds that inversion induces 
a particle-hole like transformation (MX) because levels i that 
are occupied in \xr) correspond to empty levels i in MX\xr)- 



only exceptions are states at zero energy located at a 
TRIM. These are actually the tips of Dirac nodes in the 
case at hand. Thus, Ai^^^, the difference between the 
number of even and odd states at the TRIM k,± is the 
number of ungappable states. Therefore, Au^^ cannot 
change except at a phase transition of the bulk crystal, 
where the entanglement spectrum becomes ill-defined or 
changes discontinuously. At a second order transition, 
in particular, this topological invariance breaks down be- 
cause the entanglement Hamiltonian, Eqn. (|6|, has long- 
range hopping and infinitely many bands (including the 
bulk bands) collapse to zero energy (in analogy with the 
observations of Ref . [25]) . The states may then all mix to- 
gether. We will assume here that the ungappable states 
all come in Dirac pairs^ so Ai^^^ is even. 
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The parity argument shows that there is a distinct 
phase of inversion symmetric insulators defined by having 
a protected Dirac dispersion in the entanglement spec- 
trum. What physical property distinguishes this phase? 
Clearly, surface states are not the answer since all phys- 
ical surface states are gapped. At least one distinc- 
tion is captured by the electromagnetic response of the 
system.l^^ Consider integrating out the Fermions in the 
presence of weak external electromagnetic fields. Then 
the effective action contains a term: 



Se=ie 



e 

2^ 



/ 



dTcTxY, ■ B 



(21) 



Inversion symmetry and time reversal symmetry each im- 
ply that is quantized in units of tt. Under time rever- 
sal, the sign of this term is changed since B — B. 
However, this does not rule out Se- The reason is that 
So contributes a factor of e~^^ to the weight of a field 
configuration in the path integral. The term in square 
brackets is an integer for periodic boundary conditions 
in space and imaginary time, so 9 and 6> + 27r are physi- 
cally equivalent. Hence effects due to Sq are time reversal 
invariant as long as 6> is or tt. Note furthermore that 
if the system only possesses inversion symmetry in the 
bulk, then one can argue instead that the sign of this 
term changes because E — E under inversion. 

The nonzero value, ^ = tt, is realized in strong topo- 
logical insulators, time-reversal symmetric systems with 
an odd number of surface Dirac nodes. ^ It is therefore 
natural also to suppose that an insulator with an odd 
number of entanglement nodes, but with only inversion 
symmetry, has an electromagnetic response of ^ = tt. 
(This remains true even for inversion symmetric insula- 
tors that cannot be obtained by perturbing time-reversal 
symmetric insulators. ) Note, the entanglement spec- 
trum gives us a very simple way to predict how the in- 
sulator responds to an electromagnetic field based solely 
on ground state properties. 

The physical meaning of Sq is that applying a magnetic 
field induces a parallel polarization of charge of magni- 

2 

tude Measuring this sharply is challenging, but 

it is at least in principle a physical consequence of the en- 
tanglement nodes. The locations and numbers of nodes 
are also invariant, at least without interactions, but we 
do not know the physical consequences of these proper- 
ties. 



VII. STABILITY AGAINST INTERACTIONS 

Thus far, we have discussed topological properties of 
systems without interactions. An interesting question 
is how many of the topological distinctions remain when 
the interactions between electrons are taken into account. 
When surface states exist, one can determine whether in- 
teractions affect their properties by studying whether the 
interactions are "relevant perturbations" to the field the- 
ory of the Dirac modes. ^ Furthermore, the bulk magne- 



toelectric polarizability remains quantized even when in- 
teractions are included. But the entanglement spectrum 
remains gapless even when there are no surface states and 
captures quantum numbers not accounted for by ^: at the 
very least, we believe the values of Au^^ are conserved 
modulo 2 at each transverse TRIM. (Similar invariants 
can be constructed for each direction of the entanglement 
cut, but presumably only a few of these are independent.) 

To describe inversion symmetry in an interacting state, 
one must understand how it acts on the many body states 
appearing in the Schmidt decomposition. Inversion turns 
out to be related to an anti-unitary operator, |x^) 
A4X\x'^)' The action of this symmetry is most inter- 
esting when restricted to the even-split Schmidt states, 
where half the N particles are on each side of the divider. 
Here, it satisfies (A^X)^ = — 1 when there are an odd 
number of Dirac nodes. Thus, inversion behaves exactly 
like the Kramers transformation, and each Schmidt state 
is doubly degenerate even with interactions included. In- 
terestingly, Fermion anticommutation is a key ingredient 
in establishing this fact. It has been shown recentl}'^^ 
that inversion acts on one dimensional Haldane chains in 
a similar way. 



A. Many Body States and Inversion Symmetry 

Consider the action of inversion symmetry X on the 
many body Schmidt states 



1^, 



Km)- 



(22) 



Inversion maps a right-hand state to a left-hand one. 
Since this is a symmetry, the two must have the same 
entanglement eigenvalue E^. Thus, if each eigenvalue 
has a non-degenerate eigenstate, the inversion transfor- 
mation is simple: each state maps to its partner (up to 
a phase). Things get more interesting when degenerate 
states are present in the entanglement spectrum. 

Consider a multiplet of even-split states where 
we use the label a = 1, . . . ,(i to label the d degenerate 
states in the Schmidt basis above, Eq. (22). Its image 



under inversion is a linear combination of states on the 
left: 



N 



I* 



(23) 



6=1 



If inversion is to be viewed as a symmetry, we need to 
return to the right portion of the system, to make state- 
ments about the entanglement eigenstates of a single sub- 
systems. (See Fig. jsjD.) (Symmetries mapping systems to 
one another are not so useful-for example knowing that 
the mirror image of a left-handed molecule is a right- 
handed molecule does not imply that the orbitals associ- 
ated with one of these molecules have special symmetry.) 
Inversion can be converted to a symmetry of one part 
using the pairing of left and right states in the Schmidt 
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decomposition (22). This pairing is described by a trans- 
formation M which maps |<l>^) \^a)- This rule can 
be extended consistently to an antilinear transformation 
on the rest of the Hilbert space for the left side of the 
system. 

To see this, let us describe the partner of a general 
state Ix^) on the left in a basis-independent way: 



M\x'') = {x''\PR-^\<iG)- 



(24) 



The right-hand side is a partial inner product. It 
is not a number but a wave function for the 
right half of the system since only the degrees 
of freedom on the left half are summed over: 
M\xl) = E{af } *G({af }, {af })x^*({«f }) l{«f }> 
where {<^f},{<^f} are the variables describing the two 
halves. Note that this is antilinear in |x^) because 
of the complex conjugation. Now, \x^) = Z^a^ol^o) 
is mapped by this transformation to ^^Zi^al^a)- Each 
Schmidt state, in particular, maps to its partner. Briefly, 
Eq. ( 24 ) utilizes the ground state to connect the left and 



right halves, and the operator pR~^ is there to strip off 
the different Schmidt weights. 

Now the combination of A4 and X is a stand-in for 
inversion symmetry that returns states on the right-hand 
side of the cut back to the right-hand side; in the basis 
of Schmidt states, this transformation is represented by 
KI where K is complex conjugation and / is the matrix 
given above. 

When A4X is performed twice, the state must return 
to itself, modulo a phase, hence 



(25) 



This phase factor can be only ±1 because it has to be 
reat^^. More precisely, the operator {MX)'^ has to have 
the same value (— 1)"^ for each Schmidt state which has 
half the electrons on each side, since the wave function 
Eq. is a parity eigenstate (see Appendix [B|) . 

The case / • /* = —1 is specially interesting. Time 
reversal symmetry of spin 1 /2 + n particles has the same 
property, which can be used to prove Kramers' degen- 
eracy. This algebra has no one dimensional represen- 
tation. (If / is a c-number, the product is nonnega- 
tive.) A two dimensional representation is exemplified 

by / = ^ ^ 0^ ) ' antisymmetric matrix. Since this 

is actually the generic case, whenever the inversion oper- 
ator in a particular system obeys this algebra, all levels 
must be at least doubly degenerate. 

We now show that this is indeed how the inversion op- 
erator acts in the even-split states when there is an odd 
number of Dirac nodes, by showing that S is the number 
of Dirac nodes modulo 2. For states \^r) with ^ + A: 
Fermions on the right side, (A^X)^ can be different, but 
there will still be entanglement degeneracies if 7^ 0: the 
entanglement eigenvalues for positive and negative k al- 
ways match, since the Schmidt state \^r) is degenerate 



with (and has the opposite k). This relation (valid 

for any inversion-symmetric insulator) is not as interest- 
ing as the degeneracies among the states with the same 
value /c = 0, which is unique to topological insulators. 



B. Inversion and the Topological Insulator 

We now determine the value //* = ±1 for topological 
insulators in the noninteracting limit. Note, this also 
allows us to make statements about the interacting case, 
since the value cannot jump when small interactions are 
introduced. 

The ground state of the noninteracting system can be 
built up by creating particles in all the modes Fi (see Eq. 
([l9|). Suppose there is just one Dirac node. Then 

\'^G) = lirl + li)irl+ll) 
X n (V^^4. + \/l-P«(k±)4j 

ik_L 
Pi(k_L)>i 

X n (Vp^^L + Vl -Pi(k±)4j|0) (26) 



ik_L 
J'i(k±)<| 



r 

4, = Ee"'"""/*Jr)^w^• 



(27) 

(28) 



where the first two factors are the Dirac node states. Ex- 



panding Eq. (26) gives the Schmidt decomposition. One 



should convince oneself of the following relation between 
the many-body transformations A4X and the one body 
implementation of inversion, Xs: if \xr) is one of the 
Schmidt states, and //^ is occupied in this state, then 
X^/f is empty in MX\xr) (see Fig. [stV 

Note that the second line of Eq. (26) contains the 



states that are mostly on the right, and the third contains 
those mostly on the left. Inversion maps these states 
to one another. The two states at the Dirac node are 
inversion eigenstates, so they map to themselves. Since 
they also have the same parity (say they are even, for 
instance), ra,b ^ la,b- 

The highest weight states in the Schmidt decomposi- 
tion of the wavefunction involve acting with rj's when 
Pi > 1/2 and with /|'s when pi < 1/2, so these states are 
contained in 



Ud + ll)irl+ll)\SR)\Sr.) 



(29) 



where \Sr) and \Sl) are the filled Fermi seas (all negative 
single particle entanglement energies occupied) for the 
two sides, which are exchanged by X. 

Among the highest weight states, consider the two 
states 



\pair)=rlll\SR)\SL) 



■lirl\SR)\Si 
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with an equal number of Fermions in the two sides. While 
these are converted into each other under inversion, it 
seems possible that they could mix and split, and give 
rise to one dimensional representations of the inversion 
operation. This is where showing that the inversion ma- 
trix / satisfies //* = — 1 comes in handy. 

We can get \pair) into the form of the Schmidt decom- 
position if we define the states: 

|$f)=rt|5H) |$f)=4|5L) 
\^^) = -rl\SR) \^^,)=li\SL). 
The ground state can be written as the Schmidt sum: 

(The sign s in front is an unimportant sign due to Fermi 
statistics; see appendix [bJ) 
Inversion maps 

) ^ 1*2 ) 

M simply maps |<I>f 2) -^l^f 2)5 inversion matrix 

is /* = ^ ^ 0^ ) Note that the crucial minus sign has 

arisen because of anticommutation of Fermion operators. 

Now we see explicitly that the inversion matrix satisfies 
/•/* = —1, and hence the states remain two fold degener- 
ate. As mentioned above, this result persists for any state 
with equal numbers of Fermions on the two sides. When 
interactions are included, the even-split states all mix to- 
gether, but / can be enlarged to describe the action of 
inversion on the whole space. / • /* remains equal to —1, 
so the degeneracies survive. Unlike the usual Kramers 
degeneracy for time-reversal, this result does not require 
an odd number of Fermions. 

For a more general topological insulator, one can show 
that /•/* = (— j'^l where 5 is the number of pairs of equal- 
parity single-body states with = 0, i.e. the number of 
Dirac nodes. 



The analogous quantity in the time reversal invariant in- 
sulator is the "weak" index,, ^ ^ which can be determined 
from the surface states^ or metallic topological defects in 
the crystal. No analogous physical consequence seems 
available when time reversal is broken but inversion is 
retained. 

These distinctions can be found by looking at a sam- 
ple with a finite cross-section in the y — z direction (but 
infinite in the x-direction, perpendicular to the cut). We 
have just seen that, in a noninteracting insulator, the par- 
ities of the zero- "energy" single-body states determine 
whether there are many-body degeneracies. Only the 
Dirac modes at certain TRIMs will satisfy the boundary 
conditions, and thus the value of {M.T)'^ will count the 
number of nodes at these TRIMs. 

Say the cross-section is odd x odd, with 2Ny^lx2N 
unit cells. Assuming periodic boundary conditions, the 
allowed transverse momenta will be ( (27v^+i) ^ (27V^i) ) 
(where n^, are integers whose magnitudes are less 
than or equal to Ny^ respectively). The only one 
of these which is exactly time-reversal invariant is (0,0). 
Thus all states come in pairs related by inversion sym- 
metry except for unpaired states ^(^1 5 + ^1 5) at zero 

momentum. Therefore 6 is equal to ^Ai^(o,o) and there 
is a double degeneracy if this is odd. 

To isolate another TRIM, introduce antiperiodic 
boundary conditions along one or both of the other di- 
rections. (Antiperiodic boundary conditions along for 

example, force ky± to have the form -^^^^fy, allowing 
TT but not 0.) Then double degeneracy occurs when the 
number of modes at the new TRIM is odd. 

It is possible that the actual integer value of the Ai^'s 
is conserved also when interactions are introduced, but 
there might be a more surprising classification of in- 
teracting phases. Refs.^^ gave an example for one- 
dimensional interacting topological insulators showing 
that an integer property of topological phases can be 
changed (by multiples of 8, in fact) when interactions 
are included. 



C. Distinctions Preserved by Interactions 

Now we can argue that some distinctions among in- 
sulators survive the introduction of interactions. While 
this will include the quantized electromagnetic response 
of inversion symmetric insulators, which by virtue of 
being a response function remains well defined in the 
interacting case, other more mysterious distinctions are 
also found. Without interactions, ^Au^^ is a fixed in- 
teger for each TRIM. With interactions, we will give an 
argument that suggests that at least the number of Dirac 
nodes at each a^, ^Au^^^ is well-defined modulo 2 in an 
inversion symmetric insulator. In particular, the location 
of a single surface Dirac node in the Brillouin zone sur- 
vives the introduction of interactions, as argued below. 



D. The Parity of the Many-body Wavefunction 

Many phases (like the ordered phase of an antifer- 
romagnet) have a sharp distinction only for infinitely 
large systems, but it is possible to check what phase 
an inversion-symmetric insulator is by looking at a finite 
piece of it, with an appropriate geometry. 

For a sample with periodic boundary conditions in 
the y and z directions, a finite size in the x-direction, 
and perfect inversion symmetry between its two free 
surfaces, ^Ai^(o,o) + y is the parity of the many-body 
wave-function ?/^(ri, . . . , tat) under r — r. This can 
be seen when Ai/(o,o) = 1 by inverting Eq. (29): the 



two modes corresponding to the Dirac node just map 
to themselves under inversion symmetry while \Sr) and 
\Sl) switch places. Switching them back leads to a fac- 
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tor of (-1)^/2-1 because N/2 - 1 pairs of electrons have 
been exchanged. (See Appendix[B]for more details.) This 
parity is not any easier to measure experimentally than 
degeneracies in the entanglement spectrum; a measure- 
ment would require some sort of interference experiment 
on a macroscopic crystal^^; it is just an interesting in- 
terpretation for Ai^(o,o)- (The Ai^'s at other TRIMs are 
related to the parities of samples with other boundary 
conditions. ) 

To check what phase a sample is in if it is not inversion 
symmetric, the system does have to be very large. In this 
case, one has to use the entanglement spectrum to deter- 
mine the phase, and this works only when the surfaces 
are far enough from the cut surface that inversion is an 
accurate symmetry of the entanglement surface states. 

As an aside, this implies a surprising relationship be- 
tween the magnetoelectric susceptibility and the parity 
when inversion is a symmetry: the electron- wave function 
of a crystal with a half-integer susceptibility is odd un- 
der inversion for the following geometry: the sample must 
have an even x even cross section and periodic boundary 
conditions in the y — z direction. (For the even x even 
cross sections all the TRIMs contribute to the parity). 



VIII. CONCLUSIONS 

We have studied the entanglement spectrum of insula- 
tors with non-trivial band topology. Whenever a phys- 
ical edge or surface state is present, the entanglement 
spectrum also is characterized by protected features. Al- 
though this is purely a ground state property, we showed 
it can be formally mapped to the edge spectrum of a "flat 
band" version of the physical Hamiltonian. 

On the other hand, the converse of this is not true. 
Protected features of the entanglement spectrum may oc- 
cur in systems that do not possess physical surface states. 
This can occur when a physical surface necessarily breaks 
an underlying symmetry of the bulk solid - for exam- 
ple inversion symmetry. In these cases the entanglement 
spectrum provides a unique perspective on classifying the 
phase. We illustrated this by studying three dimensional 
insulators with inversion symmetry. In general, no sur- 
face modes occur in these systems. However, since the 
entanglement cut still preserves inversion symmetry, the 
action of inversion on the entanglement eigenmodes can 
be worked out. These were shown to lead to the protected 
entanglement spectrum. 

An advantage of defining the phase in terms of pro- 
tected properties of its entanglement spectrum is that it 
allows us to deduce properties that remain stable when 
interactions are present. For the case of inversion sym- 
metric insulators, characteristic properties beyond the 
quantized magnetoelectric polarizability appear to be 
present and are stable when the particles are interact- 
ing. The corresponding physical consequences remain to 
be identified. 
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Appendix A: Inversion Symmetry in Noninteracting 
Entanglement Spectra 



The transformation Xs //^^ 



R 



can be ex- 



pressed in terms of the correlation function C. First, we 
will show how to construct the paired functions F^k^ {x) . 
(We will omit the y and z dependence.) Let us split the 
wave function into two parts, x > and x < 0, which 
can be regarded as the top and bottom halves of state- 
vectors. Then the correlation function has four parts. 



C 



Cr Clr 

^LR 



(Al) 



Since C has 1 and as eigenvalues, (7^ = (7, giving four 
matrix equations. The equation of relevance is Clr{'^ — 
Cr) = ClClR' Given an eigenfunction //^ of Cr with 
eigenvalue Pi one can obtain an eigenvector of Cr with 
eigenvalue 1 — Pi via the transformation M 

ft{x) = [M/f ](x) = ^ 



WW). 

(A2) 

The prefactor is inserted to ensure that //" is normal- 
ized. (One can check that M is a unitary transformation, 
which can be written in matrix form M 



Next one can combine //" and //^ to give wave func- 
tions for the unpartitioned system. One can show using 
the relations between the submatrices of C that Eqn. 
(19) is an eigenstate of C with eigenvalue 1, and hence 



is occupied. Now the eigenstates of Cl and Cr are the 
entanglement modes for the left and right sides of the 
systems. Hence F pairs up all the entanglement modes 
into occupied states. (There are two exceptions: if = 
or 1, then M is not well-defined. In the first case, //^ is 
definitely unoccupied, so there is no Fi corresponding to 
it. In the second case, //^ is definitely occupied and does 
not require a partner.) 

The expression for the Schmidt weights has an intu- 
itive relationship to the F's. Since each F has its own 
electron with probability one, a term in the Schmidt de- 
composition is obtained (see Fig. ^p) when a decision is 
made about which of these electrons are to reside on the 
right, say those in states Fj with j G Ar^ and which are 
to reside on the left (those in states with j G Al)- The 
Schmidt coefficients are thus given by a Bernoulli dis- 
tribution, Hj-gAh Pj Y\jeAL )' which is equivalent to 
Eq. ([l6|. (Formally, one obtains this result by expanding 

= n^k. {y/p^yl^ + vi-Pi(k^)4j io)-) 
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Since M maps the right-half to the left half of the sys- 
tem, it can be combined with inversion symmetry, /, to 
give the transformation used in Section VI, which trans- 
forms the entanglement spectrum of the right half into 
itself, Xs = IM. Since it anticommutes with Cr — ^, 
IM changes the sign of p — ^ and hence of e^, as well 
as the sign of k^. It is also unitary. Thus, this trans- 
formation has the same properties as a two-dimensional 
"CRT" symmetry of the space parallel to the cut. Here 
C is a particle- hole symmetry, is a 180° rotation and 
T is time-reversal symmetry, although the system does 
not have those symmetries independently. The C fac- 
tor changes the sign of the energy, and all three factors 
change the sign of the momentum. The product of all 
three is unitary, like X^, because particle-hole symmetry 
is antiunitary when one considers how it acts on single- 
particle states. 

This fact is enlightening since a CPT symmetry (with 
a reflection P in place of rotation) could not be used to 
prove the masslessness of the Dirac excitations: any rel- 
ativistic equation (including the massive Dirac equation) 
is invariant under CPT symmetry. (In two dimensions, 
inversion through the origin is not P symmetry, since it 
is orientation-preserving, and equal to a 180° rotation.) 



Appendix B: Parity of Finite Systems 

The result that (XM)'^ is a function only of the topo- 
logical phase of an insulator (characterized by 6) and of 
/c, the number of excess particles in a Schmidt state, can 
be proved most easily by considering the inversion trans- 
formation of the ground state wave function of a flnite 
system. This wave function must be an inversion eigen- 
state with some parity (—1)^. 

Notice by the way an important point that we have not 
discussed much: The products appearing in the Schmidt 
decomposition, Eq. ([T]), have to be antisymmetrized: 

\aR)\aL) = A.S.{^aR{ri,r2, . . . , tn^^) (B1) 
^c.L(r^+;,+i,...,rAr)}, (B2) 

where "A.S'." means to antisymmetrize in all the vari- 
ables. This wave function is not a product wave function, 
but it is the closest thing possible for Fermions: The cor- 
relations between densities on opposite sides of the cut 



vanish, surprisingly maybe'^. 

Let us call the antisymmetrized-product wave function 
^ry- The inverse of this wave function is 



= A.S. {^aui-ru. 



The inversion of the factor corresponding to the right 
side of the system X^aR = ^aR{—'^i) is a wave 
function on the left side. So we should compare 
this expression to \X^aL)\X^aR)^ which is equal to 
A.S. {<l>c,^(-ri, . . . , -rN_j^)^^R{-rN_j^^^, . . . , -tn)}. 
This differs from in the labelling of the coordi- 

nates. The parity of the permutation is (— 1)(^~^)(^+^) 
which is equal to (—1)^"^ and so 

X\^aR)\^aL) = {-1)^ \X^ aL)\X^ aR) (B3) 

Now let us calculate the parity of the ground state of a 
noninteracting insulator by looking at a particular term 
in the Schmidt decomposition. Assume as we have been 
doing that there are 6 pairs of zero energy states ^ (r|^ + 

ij^) and ^(^^5 + ill))' Since each pair has the same par- 
ity TTi = ±1, ria^b ^ 'T^iha.h' Then all the maximum 
weight states of the Schmidt decomposition are contained 

in \MW.) = nti(4-i +^^4-i)(4 +^.4)l^i^)I^L)- 



By Eq. (B3), the parity of \Sr)\Sl) is (-1)^. The 



product of all the Dirac mode operators is even under 
inversion, so the parity of \M.W.) is also (— 1) 2 . This 
also has to be the parity (—1)^ of the entire ground state 
wave function, as was claimed in Section [BJ 

Now consider a generic term of the Schmidt decom- 
position of I^g), such as \^a)- Applying inversion to 
the Schmidt decomposition maps each term to another 
term except for a factor of ( — 1)^. This is the case 

for |$«) = \^aR)\i'aL), SO (- "'^ jl^^fl) 

is another term in the Schmidt decomposition. Since A4 
is deflned to map each left-hand Schmidt state to its 
partner in the Schmidt decomposition, we can read off 
how A4 acts: it maps \^aL) to \^aR) and XI^^r) to 



(-) 



-^X\^aL)' Applying {MX)^ to I^^r) therefore 



gives ( — 1)"^+^, for a state with k extra Fermions on the 
right side. 
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In fact, every entanglement eigenstate with a finite is a 
bound state of the surface, since the bulk dispersions are 
flat and equal to ±oo. States extend farther into the bulk 
as |e^| increases. 

We need to assume that, if the model under consideration 
is a tight-binding model, then at least one center of inver- 
sion symmetry does not coincide with a site. Otherwise, 
there is no inversion-symmetric way of dividing the sites 
into two parts. 

The relationship implies /* = e^^I~^ (note, we can take the 
inverse since the determinant of / is nonzero on account of 
the initial relationship). Now multiplying by / on the left 
gives back IP — le^^ . On the other hand, taking the 
complex conjugate of the original expression gives I P — 
le-^^ so = 0, TT 

One could send the crystal toward a beam-splitting "mir- 
ror" which inverts the crystal with probability one-half, 
and then interfere the beams. Making a mirror that can 
turn the crystal inside-out like this is even harder than 
"just" maintaining the coherence of a macroscopic crystal 
during an interference experiment! 

Whether the two sides of the system are correlated or not 
on this wave function depends on whether one is looking 
at an "Eulerian representation," (as it is called in the con- 
text of fluid mechanics) with variables that are localized 
in space (such as particle number) or a "Lagrangian rep- 
resentation," with variables that are attached to specific 
particles (their positions, for example). In the Lagrangian 
representation, the variables ri and Y2 are entangled. 



